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A strategy to generate long-range entanglement in noisy quantum networks is presented. We 
consider a cubic lattice whose bonds are partially entangled mixed states of two qubits, and where 
quantum operations can be applied perfectly at the nodes. In contrast to protocols designed for one- 
or two-dimensional regular lattices, we find that entanglement can be created between arbitrarily 
distant qubits if the fidelity of the bonds is higher than a critical value, independent of the system 
size. Therefore, we show that a constant overhead of local resources, together with connections of 
finite fidelity, is sufficient to achieve long-distance quantum communication in noisy networks. 
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I. INTRODUCTION 

Quantum networks play a major role in quantum in- 
formation processing [l| , as in distributed quantum com- 
putation or in quantum communication 0, In fact 
they naturally describe the situation where neighboring 
stations (nodes) share partially entangled states of qubits 
(noisy links). One of the main tasks of quantum informa- 
tion processing is then to design protocols that establish 
entanglement between any pair of nodes, regardless of 
their distance in the network. 

Quantum repeaters offer a first solution to this ques- 
tion: one can efficiently entangle the two extremities of a 
one-dimensional lattice of size N by iterating purification 
steps and entanglement swappings [1, Q . This strategy 
needs O(log N) qubits at each node and runs in a time 
that scales as ©(polyiV). Though being very promising, 
their realization raises some technical problems, such as 
the need for reliable quantum memories @, or the diffi- 
culty in manipulating many qubits per station. The lat- 
ter difficulty is surmounted in Q, where only a constant 
number of qubits is required at each station. Various pro- 
tocols improving the rate of long-distance quantum com- 
munication have been proposed over the past few years 
(see ^] and references therein) , but either their time scal- 
ing remains polynomial in N or they are based on rather 
complicated quantum error correcting codes 

Motivated by the discovery of powerful protocols in the 
case of two-dimensional pure-state networks [To| , another 
scheme for entanglement generation over long distance 
in noisy networks was presented in [ll|. It exploits the 
higher connectivity of the nodes to gain information on 
the errors introduced by the noisy teleportations. This 
leads to a "one-shot" protocol where elementary entan- 
gled pairs are used only once, which thus relaxes the re- 
quirement of efficient quantum memories; see [l^ for the 
latest quantum communication protocol in square lat- 
tices. However, the overhead of local resources in these 
two-dimensional systems still slightly increases with TV 
(logarithmic dependence). 

In this work, we show that entanglement generation 
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over arbitrarily long distance and using the minimum 
amount of resources (constant number of qubits per node 
and quantum operations executed in a constant time) can 
be achieved in three-dimensional lattices. For this result 
to hold, the fidelity F of the elementary links has to be 
larger than a threshold F* . We first provide an analytical 
upper bound on this value and then present a numerical 
estimate based on Monte Carlo simulations. 



II. DESCRIPTION OF THE MODEL 

We consider a cubic network that consists of ver- 
tices, each of them possessing six qubits (except the ones 
lying on the sides of the cube) , on which arbitrary quan- 
tum operations can be applied perfectly. Nearest neigh- 
bors share one partially entangled state of the form 

p^{\~ ef |$+X*+| + £(1 - e) |*+X*+| 

-fe(l-e)|ci>-X$-| (1) 

Such a state can be realized as follows: a station prepares 
locally a maximally entangled pair of qubits, and sends 
one of them to a neighbor. In the quantum channel, 
the traveling qubit undergoes random and independent 
bit-flip and phase errors with probability < e < 0.5. 
This channel describes a specific physical process, but 
the generality of p is in reality complete. In fact, we show 
in App. |2]tbat any entangled state of two qubits can be 
brought to this form by local quantum operations and 
classical communication. Finally, all classical processes 
(communication and computation) are assumed to take 
much less time than any quantum operation. 

Remark 1. Physical implementations of three- 
dimensional lattices have been proposed in the context 
of quantum information pro cessing and distributed 
quantum computation [Tsl. Il4l|. For practical reasons, 
however, it may be advantageous to realize the proposed 
construction in two dimensions, using a "slice-by-slice" 
generation similar to the techniques developed in [l5l |. 
In that case note that the time required to run the 
protocol scales linearly with N . 

Remark 2. Recently, ideas of percolation theory have 
been applied successfully to the case of mixed states of 
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FIG. 1: (Color online) Non-local control phase on two qubits 
A and B, with the help of a Bell pair |'1>^)a's'- 



rank two [l6[ . In addition to the fact that the techniques 
are very different, our study is not restricted to amplitude 
damping channels, but considers full-rank mixed states 
which are robust against any small perturbations. In fact 
our protocol still works if dependent bit-flip and phase 
errors are present in the connections. 



III. A MAPPING TO NOISY CLUSTER STATES 

It was shown in [ll| how to create and propagate a 
large Greenberger-Horne-Zeilinger (GHZ) state in a noisy 
square lattice. This state is robust against bit flips if 
their rate is not too high but is very fragile against phase 
errors. Any of them indeed destroys the coherence of 
the GHZ state. Therefore, an encoding of the qubits 
is required, which leads to a logarithmic scaling of the 
physical resources per node. Since we are looking for a 
fidelity threshold, we want to create a large state that 
has the ability to correct both bit-flip and phase errors. 
Cluster states thus arise as a natural choice. In fact they 
have been shown to possess an intrinsic capability of er- 
ror correction, so that long-range entanglement between 
two faces of an infinite noisy cubic cluster state is indeed 
possible [13] ■ Our protocol is based on this construction, 
with two radical differences, however; first, the settings 
are distinct, and second, we allow only local quantum 
operations on all the nodes. 

A cluster state, which is an instance of graph states, 
can be constructed by inserting a qubit in the state |-|-) at 
each vertex of the graph and by applying a control phase 
between all neighboring pairs |18| . In our setting we can- 
not perform these control phases since they are non-local 
quantum operations, but we can add an ancillary qubit 
and perform joint measurements at each node such that 
the resulting state is a cluster state. This method has 
been described in [l^ in the case of perfect links, which 
can be interpreted as the virtual components of a large 
valence-bond state, and has been generalized to imperfect 
connections in [T^. Nevertheless, let us describe here an 
explicit (and slightly different) construction, mainly for 
completeness sake but also for relating precisely the error 
rate in the quantum networks with the one in the noisy 
cluster state. 

At each node, we add a qubit |-f) and use the noisy 
links p to indirectly perform the control phases. Let us 
first describe how this is achieved if all connections are 
perfect, i.e. their qubits are in the state |$^). We con- 
sider two nodes of the lattice, with two qubits A and B 
in the states |a) = ap |0) -I- oi |1) and |&) = 6o |0) +6i |1), 
and a connection |$+) between two qubits A' and B' , 



see Fig. [TJ We start by applying, on the qubits of the 
first node, the measurement operators 



A - |0)a(OOUa-+ |1)a(11Ua', 
Ai = |0)a(01Ua' + |1)a(10Ua', 



(2) 



with X]i=o"^l-^i ~ -"-4' which are followed by a bit flip 
X on B' if the outcome is Ai . The resulting state on A 
and B' reads |00) -I- oi |11). We then apply the second 
measurement 



Bo= |0)b(+0|b.s+ \1)a{-1\b'b, 
Bi = \0)b{-0\b'b+ \1)a{+1\b'b, 



(3) 



followed by the matrix Z on yl if we get Bi as outcome. 
Finally, A and B are left in the (entangled) state 

\Cab) = aobo |00) + aoh |01) + a^bo |10) - a^bi |11), 

which is the result of a control phase between \a) and 
\b). Clearly, if \a) — \b) = |+), the state \Cab} is the 
cluster state on two qubits. Now, let us determine which 
errors occur if we blindly perform the very same opera- 
tions but using another Bell state. It is straightforward to 
compute the results of these operations if one uses 
1$^), or l^*^) for the connections: one gets 12 0Z \Cab}, 
Z(g) I2 \Cab), or Z ® Z \Cab), respectively. Since the ma- 
trices Z commute with the control phases, it follows that 
errors do not propagate while constructing a (noisy) clus- 
ter state Pes from the cubic quantum network. Moreover, 
because of the specific choice of coefficients in Eq. ([1]), Z 
errors appear independently at the nodes. Since a node 
of the lattice has degree six at most, and two Z errors 
cancel each other, the vertices of the resulting cluster 
state suffer an error with a probability at most equal to 
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(4) 



This expression reduces to p ~ 6e in the regime of small 
error rates. Therefore, we are exactly in the setting 
of 171, where thermal fluctuations in the cluster state 



induce independent local Z errors with rate p. 



IV. LONG-RANGE ENTANGLEMENT IN 
NOISY CLUSTER STATES 

In this section, we mainly follow the construction and 
the notation proposed in |l7| . namely, the measurement 
of the qubits of pes according to a specific pattern of 
local bases. The outcomes of the measurements are ran- 
dom, but the choice of the bases establishes some parity 
constraints on them. Any violation of these constraints 
indicates an error, and a classical processing of all col- 
lected "syndromes" allows one to reliably identify the 
typical errors. This correction works perfectly for small 
error rates, but it breaks down at pc ~ 3.3% [23|. The 
difference between the present method and that given in 
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FIG. 2: (Color online) a) Bases in which qubits in the bulk of 
the cluster state are measured, b) Slightly different measure- 
ment pattern for the faces C and TZ: the central qubit is kept 
intact, and all qubits that lie in the shaded area are measured 
in the Z basis, except the ones with coordinates (e, e, 1) or 
{e,e,N), which are measured in the X basis. 



[17| is that no non-local quantum operation is allowed. 
This obliges us to design a more elaborated error correc- 
tion, leading to a different type of long-distance entan- 
glement. In fact we are not going to create a pure and 
perfect Bell pair of logical qubits, but rather a mixture 
of two entangled physical qubits. 

A. Measurement pattern and long-distance 
quantum correlations 

Let us define a finite three-dimensional cluster state on 
the cube 

C = {u = (ui, U2, Us) : 1 < ui,U2 + 1, Us < N}, 

and select two qubits A and B centered in two oppo- 
site faces C and TZ. The coordinates of these qubits are 
and {ul,u^,N), with ul^u^ + l^ {N + l)/2. 
For a reason that will soon become clear, we consider 
lattices of size = 1 (mod 4), so that ul is odd and U2 
even. Let us also introduce two disjoint sublattices To 
and Te with double spacing, where o and e stand for odd 
and even. Their vertices are 

V{To) = {u = {o,o,o)}(lC, 
T/(r,) = {w = (e,e,e)}cC, 



and their edges are given by the sets 

E{To) = {li = (o, o, e), (o, e, o), (e, o, o)} C C, 
E{Te) = {u~ (e, e, o), (e, o, e), (o, e, e)} C C. 

We also define the planes 

4"'^ ={"=(0. "2,0)} C To, 
r|"'^ = {m (mi, e, e)} C Te, 

and denote by T_J and the planes that contain A and 
B. These planes will be used to derive the Bell correla- 
tions of the future long-distance entangled state 
(we first consider that no error occurs, and then extend 
the results to noisy cluster states). Qubits that belong 
to the vertices of To and Tg are measured in the Z basis, 
while all other qubits are measured in the X basis. There 
are, however, some exceptions in C and TZ (see Fig. [2|): 
First, the central qubit is not measured, since it will be 
part of the long-distance entangled state. Second, qubits 
with coordinates ui = u\ are measured in the Z basis 
in order to create the right quantum correlations, as ex- 
plained in the following paragraph. Finally, we measure 
in the Z basis all qubits whose first two coordinates are 
(e, o) or (o, e) and which lie in the shaded areas; these 
outcomes will be important for the error correction. 

To compute the efi'ect of the measurements on the 
quantum correlations between A and B, we use the fact 
that a perfect cluster state |C) obeys the eigenvalue 
equation \C) = \C) for all u € C, where is the 
stabilizer 

Ku^Xu W Z„, (5) 

with M{u) the neighborhood of u. If we let the prod- 
ucts of stabilizers Y\u e t* ^^"^ Y\u e t- ^'^^ '-'^ 
cluster state, we find that A and B are indeed maximally 
entangled: 

XaXb \i')AB = Ax 

ZaZb \i')AB = Az Wab, (6) 

with Xx,^z G {^Ij+l}- The eigenvalues Xx.z are cal- 
culated from the measurement outcomes x and z: 

Ax = Y\. ^" II 

Az = n n 

where n^^^ = T* \ {A,B}, = T*, n^^^ = 

yK+i) u T^"^-i), and n^^^ = 4"^+'^ U T^'-'^ U 
{iul,e,l),iule,N)}\{A,B}. 

B. Error correction 

As already mentioned, measurement outcomes are ran- 
dom but not independent. It is thus possible to assign 
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FIG. 3: (Color online) a) In red (dark disks), missing syn- 
dromes on the vertices of To, in £ and TZ, corresponding to 
the measurement pattern depicted in Fig. [2J;. Rough faces 
lie on the top and the bottom of this lattice. A new coor- 
dinate system (ti,i2,t3) is introduced for the vertices of the 
sublattice. b) Same considerations for the Z correlation: the 
missing syndromes create additional rough surfaces in Te- 



on A and _B is a mixed state of the form 

PAB = FxFz |$+>(*+| +PxFz |*+X*+| 

+ Fxpz|$">($"| +PA-Pz|*"X*"h (9) 

with Fx = 1 ~px and Fz = 1—pz- This state is known 
to be distillable, and thus useful from a quantum infor- 
mation perspective, whenever its fidelity Fab = FxFz is 
larger that one- half [23] . This can be achieved when the 
error rate p is smaller than a threshold p* . In the next 
paragraphs, we first prove a lower bound on this value, 
p* > 1.17 X 10~^, and then present numerical results, 
showing that the real threshold is indeed much larger: 
p* > 2.27%. 



1. Correlation loss due to the missing syndromes in To 



to most vertices Ui € Ti, with i = o or e, the parity 
syndrome 

s{ui) = n 2;„ z^, (8) 

where Afi{ui) designates the neighborhood of Ui in Ti. 
Since this equation arises from a product of stabilizers, 
YivGJViiui) -^v, we have that s{ui) — 1 if no error 
occurs on the qubits of J\fi{ui). The key point of the 
construction is that a Z error on any edge of Ti changes 
the sign of the two syndromes at its extremities. This 
is due to the fact that Z errors do not commute with 
X measurements, while outcomes z are not affected 
by them. The sublattices are treated separately, but 
in a similar way. We refer the reader to [ItL [21| for a 
detailed discussion of the error recovery or to App. |B] 
for the basics to understand our protocol. In contrast 
with [l7j . and apart from the rough faces present in 
any surface code, we also suffer a lack of syndrome 
information in C and TZ. We cannot have a perfect and 
complete syndrome pattern for both To and T^ in these 
faces; for this to happen one should be able to measure 
both X and z eigenvalues of the concerned qubits, which 
is impossible, or apply non-local quantum operations, 
which we do not allow. Actually, useful long-distance 
quantum correlations can still be created if one performs 
the measurements depicted in Fig. [5]d: half outcomes 
are used to gain information on To, and symmetrically 
for Te, see Figs. Hand IH 

As an example of the effect of the unknown syndromes 
in C, let us consider that an error occurred on the center 
qubit A, and that all other qubits did not suffer any 
error. Since we do not know the syndromes of To that lie 
directly below and above TJ, we are not able to restore 
the X correlation. This occurs with probability px — 
p + 0{p'^). From this fact, one finds that the final state 



Paths of errors, which we generically denote by F, have 
a non-trivial effect on the X correlation if they cross the 
plane an odd number of times, as depicted in Fig. |4l 
Moreover, the number I of errors which actually occur on 
a path F is at least L/2, where L denotes its length. This 
is the case because our error correction always leads to 
a minimum pairing of the syndromes s — ~1. We now 
follow Chap. V in [2l| to find an upper bound on the 
probability px of inferring the wrong quantum correla- 
tion: 

< 2 ^ prob(F££) + prob(F£7j) 

Tec r^Tz 

+ ^prob(Fre)+4^prob(Fr£), (10) 

where B and T stand for the bottom and top faces. Note 
that we already took into account the symmetries of the 
problem in this expression. For convenience, let us now 
introduce a new coordinate system {ti, ^2,^3) for the ver- 
tices of To, such that -No < h < No, -No <t2 < No, 
and < ^3 < 2No, with No ^ {N - l)/4, see also Fig. [1 
In this coordinate system, paths of errors Ft-^ travel a 
distance L > No and can start from No missing syn- 
dromes in To (lower triangle in C). Because for each 
vertex there are, in a cubic lattice, at most 5^/2 self- 
avoiding walks pointing upward, we find that the last 
term of Eq. (IIOD is upper bounded by 

E p-b(rr.) < E Y E 

Ftc L>N^ l=\L/2] ^ ' 

where [i/2] denotes the smallest integer not less than 
L/2. The sum over together with the binomial coeffi- 
cients, counts all possible paths of errors that appear in 
a given walk. One can check that the bound tends to 
in the limit No^ 10^p(l -p) < 1, i.e. \i p < 1%. 
The same result holds for the paths F^tj and Ft-^; note 
that this value is about three times smaller than the real 
critical point po ~ 3.3%. Similar considerations for the 
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FIG. 4: (Color online) Some paths of errors that have a non- 
trivial effect on the long-distance entanglement: any path 
stretching from one shaded area to another, and crossing the 
plane an odd number of times, degrades the X correlation 
between A and B. The shaded areas, which partially wrap the 
cube, are of two types: the top and bottom ones are the usual 
rough surfaces present in (three-dimensional) surface codes. 
The left and right shaded areas represent the unknown syn- 
dromes of To, see Fig.|3^. The situation for is very similar 
(the picture is rotated by 90° around the AB axis), with the 
difference that all shaded areas are rough surfaces in this case. 



Fab 




FIG. 5: Fidelity Fab ~ FxFz of the long-distance entangled 
state pAB as a function of the error rate p (in percent). The 
critical value F* = 0.5 is reached at p* = 2.27 ±0.03 [%]. The 
upper curve represents the probability of success Fz of the 
error correction in Te, while the middle one is the function 
Fx- Corresponding series expansions for p ^ 1 are plotted 
with dashed lines. 



paths finally yield, for p < 1%, 
Px < 2^prob(r cc) 

t2>l L>t2 l>[L/2] ^ ^ 

This bound never tends to zero, but still converges if p 
is small enough. Before computing a threshold for Fab, 
however, we first have to consider the errors made in the 
other sublattice. 



2,. Loss of correlation in Te and fidelity of the final state 

The situation for the Z correlation is very similar to the 
previous case, since the measurement pattern is symmet- 
ric. Nonetheless there is a small difference: the missing 
syndromes do not lie on the vertices of the sublattice, but 
rather on its outer edges C and TZ. This creates additional 
parts of rough faces in Tg. It follows that the correspond- 
ing paths of errors have their first and last edges pointing 
in the direction, so that a slightly better bound for the 
error pz can be derived: 

pz<2E2iiE^E [''^p'{l-pf-'■m 

ti>l L>ti+2 l>\L/2'\ ^ ^ 

Combining the two bounds on px and pz we find 

p< 1.17X 10"^ i^AS > 1/2, (13) 

which corresponds, via Eq. to an error rate e < 
1.95 X 10~* in the initial connections. This value is quite 
small, mainly because our counting of paths of errors is 
very crude. Note that there are only few such paths of 



small length, and therefore this analytical bound could 
be increased by carefully computing its smallest orders 
in p. At this point, however, we prefer to turn to Monte 
Carlo simulations to find a much better estimate of the 
error threshold. We refer the reader to App. [C]for a de- 
scription of the algorithm; in particular we propose an 
intuitive and efficient method, even if not optimal, to in- 
fer the value of the missing syndromes. The result of 
these simulations is plotted in Fig. [5l long-distance en- 
tanglement is achieved for error rates smaller than 

p* « 2.27%, (14) 

i.e. e* ~ 3.86 x 10"'^ for the original lattice. Before con- 
cluding, let us comment on these thresholds: 

• The values of the unknown syndromes are not opti- 
mally inferred in our algorithm, and therefore a higher 
value of p* may be found. However, it is clear that it 
cannot exceed the critical error rate Pc ~ 3.3%. 

• One could get a higher threshold e* by directly com- 
puting Fab as a function of e. In fact, errors in the 
faces £ and TZ do not appear with probability p 6e, 
but only with probability p ^ be. 

• Our measurement pattern puts To and Tg on the same 
footing (Figs. [5] and [3]), but it could be profitable to 
get more information on the unknown syndromes of To 
since the X correlation is more sensitive to errors. 

• Finally, as suggested in [13, Rem. 2], lattices of size 
\og{N) X log(iV) X A'' may also be appropriate for gen- 
erating long-distance entanglement. This result also 
holds in our setting, because additional errors only ap- 
pear in the faces C and TZ and not in the bulk of the 
lattice. 
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V. CONCLUSION 

We have investigated the problem of generating 
long-distance entanglement in noisy quantum networks. 
We have focused on three-dimensional regular lattices, 
whose edges are full-rank mixed states of two qubits. 
We have proven that entanglement can be established 
between two infinitely distant qubits if the fidelity of 
the connections is large enough. Our protocol starts 
by transforming the quantum network into a thermal 
cluster state. Then, all but two distant qubits are 
measured according to a specific pattern of local bases, 
and a syndrome-based error correction is performed. 
The error recovery is very similar to the one used for 
planar codes, with the difference being that our setting 
does not allow one to get complete information on the 
syndromes. Nevertheless, useful quantum correlations 
can be created between the two unmeasured qubits if 
the error rate is smaller than a critical value. We have 
given both an analytical lower bound on this value and a 
numerical estimation (about 2%) based on Monte Carlo 
simulations. 

In conclusion, we have shown that a constant over- 
head of local resources is sufficient to achieve long- 
distance communication in quantum networks. This con- 
trasts with previous one- or two-dimensional strategies, 
in which the physical resources per station increase with 
the distance. Our protocol requires perfect local quan- 
tum operations, which is somehow justified by the fact 
that most errors occur while sending quantum informa- 
tion between stations. It would nevertheless be of fun- 
damental interest to design fault-tolerant protocols, in 
two or three dimensions, for which the overhead of local 
resources is as small as possible. 
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Appendix A: Elementary entangled pairs of qubits 

We show in this appendix that there is no loss of gener- 
ality in choosing the elementary links to be described by 
Eq. ([1]). First, it is well known that any two-qubit entan- 
gled state can be brought to the rotationally symmetric 
mixture 

Wf^f |$+X*+| + l*+X*+l 

+ ^^I*-X*-| +^^I*-X*-|- (A.i) 

In fact, this is achieved by applying random bilateral 
rotations, locally, to each qubit. This state is called a 



Werner state [2^ and has the same fidelity F as the initial 
state from which it derives. Let us denote Wp by its 
components in the Bell basis, Wp = {F, ^~W~^' 
and define the two unitaries 

One can check that the separable operation H (g) H 
applied on a Bell-diagonal state switches its coefRcients 
and 4'"'", while the coefficients and ^~ are 
unaltered. A similar result holds for H' (S) H', which 
only switches the components and '^'^ . Suppose now 
that an entangled pair Wp, with F ^ 1 — 3e^, has been 
created between two neighboring nodes. This already 
sets the coefficient ^E"" to the desired value e^. Then, 
apply H' (g) H' with probability p = j^^, and 1 (g) 1 
with probability 1 — p on Wp ■ This leads to the state 
((1 — e)^,e(2 — 3e),e^,£^). Finally, repeat the operation 
by applying H ® H with probability one-half. Both 
coefficients '^'^ and are set to e(l — e), which proves 
that the state p given in Eq. ([T]) is indeed general. 



Appendix B: Basics of syndrome-based error 
correction 

Let us consider the sublattices To and Tg described in 
the text, in which Z errors occur independently on each 
edge with probability p, and assign to each vertex the 
syndrome s = if it is connected to an even number 
of erroneous edges, and s = — 1 otherwise. In the case of 
perfect and complete syndrome information, one knows 
exactly where all paths of errors start and end: this oc- 
curs at syndromes s = — 1. In the regime of small error 
rate p, it turns out that the best error recovery strategy 
is to pair these syndromes such that the total length of all 
pairings is minimized. Then, one connects any two paired 
syndromes by a path of minimum length, and artificially 
introduces Z "errors" along these paths. This creates 
loops of errors in the cluster state, which, however, do 
not cause any damage to the long-distance quantum cor- 
relations. In fact these loops cither do not intersect the 
planes and or cross them twice, and consequently 
do not modify the eigenvalues in Eq. (jH]). 

Problems arise because some syndromes are unknown. 
For instance, consider the edges that have only one ex- 
tremity in V{To) or V{Te): their coordinates are (o, 0,o) 
and {o,N — l,o) in To, and (l,e, e) and {N,e,e) in Tg, 
see Fig. [31 These are the rough faces described in [13 , 
and errors on these edges change the sign of only one 
syndrome (and not two) in the corresponding sublattice. 
An equivalent viewpoint is that both extremities of these 
edges indeed belong to To or Te, but we do not have 
access to their outer syndrome. The consequence of this 
lack of information is that some paths of errors F are 
not closed anymore, but rather originate from a missing 
syndrome and terminate at another. Typically, these 
open paths enter only superficially the lattice if the error 
rate p is small, but they start stretching from one side 
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FIG. 6: (Color online) An example of how unknown syn- 
dromes, which here lie on the boundary, are assigned the value 
±1. We first find all odd-size clusters d of syndromes — 1, 
which are drawn in light yellow (light gray), and pair them 
by increasing distance. A cluster may be left alone; in that 
case we add it in the end of the sorted list of cluster pairs Pi. 
Then, for each Pi, we check if it is favorable to add two new 
syndromes —1. Here, this is the case only for P4 and P5. 



to another as soon as p exceeds the value Pc ~ 3.3%. In 
the latter case, paths of errors can cross an odd number 
of times the planes of correlations, which results in a 
complete loss of long-distance entanglement in the limit 
N ^ 00. 



Appendix C: Monte Carlo simulations 

We now describe the algorithm used for computing the 
data of Fig. [SJ The correction procedure is very similar 
for the two sublattices To and Te, and consists of two 
main parts. First, given a lattice with random errors, 
we infer the value of the syndromes for which we have 
no information. Second, we proceed with the usual er- 
ror recovery. The program outputs 1 if the correction is 
successful, and otherwise. 



a. Inferring the missing syndromes 

We propose a very simple way of assigning the value 
+1 or —1 to the missing syndromes, so that a good ap- 
proximation of the optimal configuration is found. To 
that end, it is helpful to consider a typical realization of 
a noisy cluster state in the regime of small error rates, as 
depicted in Fig. [51 Our algorithm reads 

• Initialize all unknown syndromes to +1. 

• Using nearest-neighbor site percolation, find all clus- 
ters Ci of syndromes —1. Keep only the clusters of odd 
size, and for each compute the minimum distance di to 
a closest unknown syndrome Si. [1^ Let n denote the 
number of such clusters. Note that we do not consider 
clusters of even size, since good pairings can be found 
for them, individually. 



• Calculate the distance dij between all pairs of clusters 
Ci and Cj. This distance could be the length of the 
shortest path from Ci to Cj , but in practice it is much 
easier to calculate the distance between their "centers 
of mass." 

• Find Ca and Cb such that dab = min{dy }, and create 
a pair Vi = {Ca, Cb}. Remove Ca and Cb from the list 
of clusters, and repeat the procedure until no cluster 
is left. In case of odd n, add an extra "pair" {C,, C^} 
for the remain cluster, with da — od. This creates the 

list {Pi,...,P„/2}. 

• For each Pk = {Ci,Cj}, check if dij > di + dj. If this 
inequality holds, inverse the value of the corresponding 
missing syndromes: Si (—Si) and Sj 4— {—Sj). 

The proposed algorithm is optimal in the regime of 
very dilute errors, but this is not true for high error 
rates anymore (even if results are good for all p). Note 
that there exist optimal algorithms which are based 
on minimal perfect matchings in weighted graphs and 
are run in a polynomial time (see [2J|, Chap. 4). For 
three-dimensional lattices, however, the number of 
edges in these graphs scales as 0{N^) [they are nearly 
complete graphs on 0{N^) vertices], and therefore these 
algorithms are not so efficient in practice. Nevertheless, 
it would be very interesting to implement an optimal 
algorithm and decide whether the unknown syndromes 
are responsible for the threshold, or whether the equality 
p* = Pc holds. 



b. Error recovery 

We use the well-known and efficient algorithm de- 
scribed by J. Edmonds in [25l | to find an optimal pairing 
of the syndromes —1. The error correction is successful 
if the parity of paths of errors crossing the plane of 
correlation is even. Simulations of the error corrections 
have been performed for various lattice sizes (up to 
sa 15^ nodes), and for both X and Z correlations. The 
extrapolation to infinite lattices is done by fitting the 
data with an exponential function, see Fig. [T] Results 
are plotted in Fig. [S] 

Finally, let us present evidence that our algorithm gives 
correct and optimal results in the regime of small error 
rates. Considering the series expansions of px at first or- 
der in p, one sees that only three edges of £ may degrade 
the X correlation: these are the bonds in To that cross 
and whose first coordinate ti belongs to {—1,0,1}, 
see Fig. [3^. At second order, one can check that the 
probability to infer the wrong X correlation due to the 
missing syndromes in C is pj^ — ip + ASp^. Therefore, by 
symmetry, the fidelity Fx — 1 — Px{l — p^) — p^{l — Px) 
reads 

Fx = l-6p~ 78p2 + 0(p3). (C.l) 
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It is easy to see that a single error in Tg cannot damage 
the Z correlation, and a careful counting of configurations 
with two errors yields — 19p^. Consequently we find 



0.70 



0.65 



No 



Fz = I-38p2 + 0(p3). 



(C.2) 



FIG. 7: Fidelity Fx in the limit No oo, for a fixed error 
rate p. We use the function + ae^*'^° to fit the data 

computed from lattices consisting of (2A^o + f ) x {2No) x {2No + 
I) nodes. At least W simulations have been run for each 
value of No £ {2, 3, . . . , 7}. In this example the error rate is 
p = 2.2%, and the fit yields = 0.699 ± 0.003. 



These two series expansions are plotted (dashed lines) 
in Fig. [5j they agree perfectly with the results of the 
Monte Carlo simulations and thus validate our algorithm. 
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